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Abstract 

In analogy with an experimental setup used in liquid helium, we use a pointlike probe to study 
superfluids which have a gravity dual. In the gravity description, the probe is represented by a 
hanging string. We demonstrate that there is a critical velocity below which the probe particle 
feels neither drag nor stochastic forces. Above this critical velocity, there is power-law scaling for 
the drag force, and the stochastic forces are characterized by a finite, velocity-dependent temper- 
ature. This temperature participates in two simple and general relations between the drag force 
and stochastic forces. The formula we derive for the critical velocity indicates that the low-energy 
excitations are massless, and they demonstrate the power of stringy methods in describing strongly 
coupled superfluids. 
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1 Introduction 



Liquid Helium, at temperatures below ~ 2.17° K, is the paradigmatic example of a superfluid: a 
liquid with a broken global symmetry. Since its Helium atom constituents are strongly interacting 
it is difficult to obtain good theoretical control over the properties of the superfluid beyond the hy- 
drodynamic approximation. This is exactly what the gauge-string duality provides: recent works 
ni 121 El IS [51 m offer theoretically tractable examples of a strongly coupled superfluid, holograph- 
ically dual to superconducting black holes in five-dimensional anti-de Sitter space (AdSs) via the 
gauge-string duality |I7][8l|9l|. 

In this work we will use a holographic model of a superfluid to study its response to a heavy, 
pointlike probe. Such an arrangement is similar to an experimental setup where a heavy ion is 
dragged through superfluid helium [[TOll . Due to the low viscosity of the superfluid the ion moves 
without friction until it reaches the critical velocity for quasi-particle creation upon which it slows 
down. In superfluid helium, the quasi-particle excitations are rotons or vortices, depending on the 
details of the experimental setup. We will show that in the holographic superfluids we study, there 
is also a critical velocity which can be determined from the bulk geometry but now the excitations 
associated with the onset of drag can be thought of as massless. We also study the stochastic force 
acting on the moving probe and extract from it the transverse and longitudinal mean momentum 
transferlll 

Holographic superfluids are different from liquid helium, being based on strongly coupled 
gauge theories. Moreover, the infrared dynamics of the superfluids we study is characterized by 
emergent conformal symmetry. The critical velocity that we compute from the low-temperature 
thermodynamics is the one that the emergent conformal group leaves invariant. It defines an index 
of refraction of the holographic superfluid. Thus we encounter a surprising confluence of ideas, 
weaving together superfluidity, emergent symmetry, black hole physics, string dynamics, and re- 
fractive signal propagation. 

2 The background geometry 

The gauge-string duality allows one to map a theory of gravity in an asymptotically five dimen- 
sional AdS space to a strongly coupled, large A^, SU (N) conformal field theory. Our starting point 

'The formulas for the mean momentum transfer have been obtained independently in ifTTI in a different context, 
and a significant simplification of the expression for kl was carried out using techniques borrowed from that work. 
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is the five dimensional bulk action S = ^ J d^x^/^ C with 

^ = ^ - -/l^ - \ {{Q.nf + {d,d - qA,f\ - V{v) • (1) 

From here on, we will set = 1/2. In [|51 it was shown that ([T]) can be constructed from a 
truncation of HE supergravity on a Sasaki-Einstein manifold. (For a related example in M-theory, 
see [[El [131.) In this truncation, one finds that ^(r/) = sinh^ rj, V{ri) = — cosh^ ^{5 — coshr/), 
q = Vs/L and the real scalars rj and 9 describe a complex scalar whose target space is the Poincare 
disk. This complex scalar is charged under the f/(l) gauge field (where F = dA), so non-zero 
r] breaks the gauge symmetry. If a solution to the classical equations of motion of © is found, 
then it can be lifted to a solution of type IIB supergravity. As explained in [[H [2|| the breaking of 
the U (1) gauge- symmetry in the bulk corresponds to a formation of a condensate in the boundary 
theory. 

We are interested in the bulk solution which corresponds to a zero temperature configuration 
where the charged scalar has no non-condensed component. In lfT4ll it was suggested that the zero- 
temperature limit of black hole solutions of theories similar to O would be a domain wall with 
anti-de Sitter space on both sides: If we parameterize the line element, gauge field, and scalar 
field as dslj = e^^^'') {-h{r)dt^ + dx^) + dr^/h{r), A^dx^ = ^{r)dt, i] = r]{r), 9 = then the 
blackening function h has positive slope everywhere, is approximately 1 in the ultraviolet (UV) at 
large positive r and asymptotes to a constant, which we denote Vj^, in the infrared (IR) where r is 
very negative. The warp factor A{r) also has positive slope everywhere and interpolates between 
two different linear functions of r, namely r / Ljr in the IR and r/L in the UV. Similarly, the scalar 
r]{r) rolls down from a local maximum of the potential V in the ultraviolet to a local minimum 
in the infrared. Note that Lir and rji^ are related through \^(?7ir) = We will interpret 

this domain wall solution, in the dual gauge theory, as the ground state for a finite density of [/(I) 
charge, in which all the charge is carried by a superfluid condensate. 

In the case of the IIB Lagrangian discussed above, ?7ir = cosh^^ 2 and rjuy = 0. An explicit 
domain wall solution for this theory was constructed in [[6l, where it was found that fjR ~ 0.373. 
This domain wall is a typical example of a zero temperature configuration with conformal invari- 
ance in the infrared as well as the ultraviolet. The result of [6j has been reproduced in figure [T] 
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Figure 1 : A plot of h{r) for the particular type IIB theory described in the text, at zero temperature. 
The scalar flows from an AdS^ region in the UV to an AdSs region in the IR. 



3 The drag force 

With the bulk description of a zero temperature superfluid in hand, we examine a heavy pointlike 
particle moving through the superfluid at constant velocity. On the gravity side this is described 
by a trailing string, similar to the treatments of IfTSl [T6ll . One end of the string is located at the 
asymptotically AdSs boundary and specifies the location of the probe particle in the boundary 
theory. The other end of the string trails off into the deep IR. The momentum flowing down 
the string is equal in magnitude to the drag force acting on the boundary theory particle. Thus, 
to find the drag force one needs to look for an appropriate string configuration and compute the 
momentum flowing down the string. 

We will work with a string action of the form 

Sstring = j dadr Cstring whcrc ^string = ~ 2^ ^ ^^^^ ~ '^^^ ^c,X^'^pX''g^y . (2) 

Here g^y is the five-dimensional Einstein frame metric and (r, a) are worldsheet coordinates. In 
the particular case of the IIB action considered in [[51 one finds that Q = cosh |. The procedure of 
extracting the drag force from the action is straightforward and here we merely quote the result: 

^drag = ^Ve{v - Wir) , (3) 
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Figure 2: A plot of the trailing string with one endpoint on the boundary and the other in the IR 
region of AdS space for the IIB theory of [5]. For v < vm there is no solution to /i(r^,) = v"^ and 
the only possible configuration is a string hanging straight down and moving at constant velocity. 



where A* = A{r^), = Qijjir^)), and so forth, and r^, is determined implicitly by the equation 
h{r^) = . The result © (which was basically already understood in [fTTl ) provides some inter- 
esting intuition about what warped geometry means: the string frame warp factor on the spatial 
dimensions x^, x^, and is e^^Q, and, up to velocity-independent constants, this string frame 
warp factor evaluated at r = r^, (when r* exists) is the coefficient of kinetic friction for a heavy 
pointlike probe. The surprising feature of ^ is the appearance of the step function 9 which van- 
ishes when its argument is negative: Once the velocity of the probe is smaller than v\yi then = 
has no solution, and the drag force vanishes. This is depicted in figure [U By Taylor expanding the 
drag force, the asymptotic behaviors of Fdrag as v approaches 1 from below or f ir from above can 
be written in terms of the infrared or ultraviolet parameters of the theory. 



-^drag 



where (Too) is the energy density of the zero temperature configuration in the boundary theory, 
A$ is the dimension of the U{1) current in the infrared defined so that a conserved current has 
dimension 3, and bh is an extra parameter which characterizes the deviation of h from a constant 
value in the IR. For the type IIB theory of f5l, A$ = 5. The details of the derivation of © and Q 
can be found in appendix 1X1 
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The onset of the drag force in the holographic setup we've just described seems to be associated 
to a phenomenon analogous to roton emission and to Cherenkov radiation — the latter providing a 
slightly closer analogy, as we now explain. Low-energy photons in a refractive medium comprise 
a free conformal field theory with a reduced speed of propagation relative to the speed of light in 
vacuum. When an electrically charged particle travels through the medium faster than the speed 
of light, nearly free photons are emitted in a process called Cherenkov radiation. What probably 
happens in our setup is that an operator which couples to the pointlike probe mediates strongly 
interacting "radiation" into the medium, and this radiation propagates at a velocity ^^irE Although 
this line of thought is plausible, it surprises us that the scaling F^rag ^ {v — wir)" for v slightly 
larger than t>iR is characterized by an exponent a = (A$ — 4)~^: The spectral measure of the two- 
point function of the operator that the probe couples to is the analog of phase space for photons 
emitted through the Cherenkov process. Naive phase space measure arguments applied to the case 
of weakly coupled quasi-particle emission suggest instead a scaling exponent a ~ A, where A is 
the dimension of the operator responsible for the quasi-particle emission. An extension of these 
phase space arguments to strongly coupled theories, a la [|T8l, still gives an a ~ A scaling where 
A would now be the dimension of the "quasi-unparticles." Perhaps a better understanding of how 
ultraviolet cutoffs affect the strongly coupled radiation rate would resolve the quandary of the 
scaling of a; or perhaps there is an operator mediating the radiation whose anomalous dimension 
in the infrared is inversely proportional to A$. This issue certainly warrants further attention. 



4 Thermodynamics and the critical velocity 

In [[T9ll it was shown that a relativistic gas of particles in the grand canonical ensemble with dis- 
persion relation u = kcq leads to a relation of the form = sT / {sT + fxp) where s is the entropy 
density, T is the temperature, /i is the chemical potential, and p is the charge density of the quasi- 
particles. We now show, using arguments entirely on the gravity side, that a similar relation holds 
in the low-temperature limit for solutions of Qi: 

Vm = lim -7^^ , (5) 

where p„ is the density of charge in the normal state, defined as in EOl by the thermodynamic 
relation dP = s dT + p„ dp — ^d^'^, where = and (p is the Goldstone boson of the broken 
U{1) symmetry. In our derivation, we will not require the assumption that there are low-energy 

^We thank D. Son for discussions on this point. 
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quasi-particles that propagate with speed wir. 

To prove ([5]) we use the formalism developed in ||20l . So far we considered the zero temperature 
solution of the theory. If the temperature is finite then the solution will be a black hole, by which 
we mean that the infrared limit of the theory is characterized in the bulk by a black hole horizon 
at r = To where /^(ro) = 0. The transition from a black hole solution to a domain wall solution 
occurs by the blackening function h developing a "double shelf" structure. That is, the horizon 
is pushed to the deep IR, before which h is characterized by a region where it is approximately 
constant at a value approaching t>jj^. The double shelf structure was first discussed in [211 • 

In these positive but low temperature solutions, we consider vector perturbations of the station- 
ary metric of the form A,j.i = a(r) and g^^^ = e'^^^'^^g{r), and solve the linearized equations of 
motion for g and a. The explicit form for the equation of motion for g can be written as a total 
derivative by making use of the equations of motion for the background solution. After integrating 
it once it takes the form e^'^g' + e^^a^' = Q2 where Q2 is an integration constant. By expanding 
the fields a and g near the boundary and using the relativistic two-fluid model for superfluids, intro- 
duced in |22j and elaborated upon in [l4l|23l[20l, one can show that Q2 = sTu^^ + {u^^ — Vx^)^Lpn- 
Where u and v are the velocity of the normal component and the superfluid velocity respectively. 
This leads to Q2 = sTu^i + {u^i — t^^O/^Pn- Making use of Lorentz invariance, a similar analysis 
can be carried out for the equation of motion for a. One finds he'^"^{a(f)' — a'cj)) = hQ2 — gQi + Qz 
with Q3 = {vx — Ux)fipn- To obtain ([5]) we can evaluate the equation of motion for a on the "in- 
frared shelf" described earlier, and then take the zero temperature limit. The finer details of this 
analysis are left to appendix |Bl 

Considering ([5]) together with the computation of |fT9l and the expression for the drag force 
for V near vib., it is tempting to speculate that the low temperature thermal excitations are massless 
particles. However, one should keep in mind that our system is strongly coupled, so it's more likely 
that the critical velocity is associated with a strongly coupled "soup" of excitations whose collective 
behavior is in some ways similar to that of massless particles, but cannot be fully described in terms 
of quasi-particles. 

5 Stochastic forces 

Having determined the drag force acting on the probe particle and related the threshold velocity vi^ 
at which it becomes non-zero to thermodynamic quantities, we proceed to compute the stochastic 
forces on the probe. In a Langevin formalism with white noise, one expresses the total force acting 
on the probe as dp/dt = Ftot = -^drag + -^l(^) + FT{t), where FL{t) is the stochastic force in the 
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same direction as the motion of the probe particle while -Ft(^) is orthogonal to it. The expectation 
values of Fl and Ft vanish, while (ts)) = t^Lp'p'^iti - ^2) and {F^{ti)F^{t2)) = 

KxiS^-' — p'p')5{ti — t2), where is the unit vector in the direction of the momentum p. 

Li Il24ll25ll26l it was shown explicitly that the stochastic force F^ + Ft owes to the existence of 
a horizon on the string worldsheet: the endpoint of the string on the boundary is causally accessible 
to only that part of the string worldsheet with r > r*, where is determined as before from = 
v"^. (If V < fiR, then there is no worldsheet horizon. We will focus on the case wir < v < 1 from 
now on.) The metric on the string worldsheet above the horizon can be written in diagonal form 
ds^g = ■jabdcr^'da^ = —e^^{h — v'^)dT'^ + + ^T-S^ ) ^''^^^ where the worldsheet coordinates 
(r, r) are specified by the indices a, b. Note that r = t only on the boundary. The worldsheet 
Hawking temperature is given by 



Following arguments similar to the ones leading to ([4]), one finds 



(6) 



Tws 



ML, 



V 1 



(7) 



provided Q approaches a constant in the UV and IR. When A<^ = 5 then Tws is linear in t> — t>iR 
for V close to t^iR. 

To study the stochastic forces further, one should consider perturbations of the string around 
the trailing string solution which lead to the drag force result ([3]). Computations of this sort where 
carried out in [|25]|26ll27ll28l and can be extended to general holographic backgrounds of the form 
we've discussed. The technique for carrying out this extension is to use the method of matched 
asymptotic expansions for approximate expressions for the fluctuations at small frequency and in 
the deep infrared. The result is 

kt = 



(9 log V 



The details of this computation are left to appendix O We note that the form of is interestingly 
similar to the Einstein relation /£|™stcin _ —2Fdra.gT/v, which can be derived by demanding that 
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stochastic evolution based on the total drag force Ftot equilibrates to a thermal distribution of 
temperature T, see for example ||29l . But the differences between kt and 

^Einstein gj.g 

significant: 

besides constraining kt rather than k^, the involvement of T^s in the expression for kt makes it 
hard for us to see what the late-time equilibrium distribution of momenta would be when the probe 
mass is made finite. 
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A Details of the drag force computation 



Our starting point is the string action given in Consider the following ansatz for the string 
embedding. 

/ r + C(r) \ 
VT + v({r) + C,{r) 





v 



(9) 



where we have set a = r. In the original treatments [[T5l[T6l . the function ((r) was not introduced; 
indeed, the drag force computation is straightforward with ((r) set to 0. But, as pointed out in 
OOH . the choice 

CM--/ (10) 

makes the worldsheet metric diagonal, which in turn makes the treatment of fluctuations easier. In 
(flOl ) and hereafter, primes denote d/dr. For the sake of consistency, we will work throughout with 
C(r) as given in (flOl ). 

Since ^' and not ^ enters into £string, 



(11) 



is constant when the equations of motion of the string are obeyed. An analysis of the worldsheet 
current of spacetime energy-momentum analogous to the one in |fT6l shows that is equal to the 
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drag force -F^rag- Solving (fTTI) for ^'(r) leads to 



,'(r\ = -lL- (12) 
' he^y he^^Qy{27ra'y -71^ ^ ^ 



Requiring that the solution is real everywhere and that only one string endpoint lies on the AdS 
boundary, we are lead to either tt^ = 0, or 



27ra' 

where is determined implicitly by the equation 

h{n)=v\ (14) 

As in the main text, we will henceforth use a subscripted star to denote functions evaluated at 
r = r*. An important subtlety in (fT3l) is that if v < fiR, then there is no solution to (fT4l) . In that 
case, the only solution to the classical equation of motion that extends arbitrarily far down into the 
bulk is the configuration with ^ = ^ = 0. If f > t;iR then the vr^ = solution is no longer physical 
because then the worldsheet metric is not Lorentzian for all r. Putting everything together, we find 
the drag force dS]). 

The asymptotic behavior of -Fdrag as v approaches 1 is easy to understand based on the fact that 
the leading departure of the background geometry from pure AdS 5 is 

h{r) = 1 - ^e-^^« + . . . , Air) = ^ + . . . (15) 

where (Tqo) is the energy density of the dual field theory configuration. Using this expansion in 
© leads to ©. Recall that we have set Kg = 1/2. 

The behavior of Fdrag as v approaches vi^i depends on some detailed properties of the infrared 
geometry. The key parameters are the infrared dimensions and A$ of the operators dual to 
7] and Afj_; the bulk scalar and gauge field obtain an effective mass which can be computed by 
linearizing the equations of motion around the infrared fixed point. For the particular geometry in 

A^ = 2 + 2v^ A$ = 5. (16) 
For general values of A^ and A$ one finds the following solutions to the linearized equations in 



"^In 1141 [3TI . A$ was chosen differently. Our conventions here are such that A$ = 3 when the photon is massless, 
which is when it corresponds to a conserved current in the dual field theory. 
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the IR: 

h = v^^{l + 6,,e^'^"/^™ + . . .) 
V = ViR + se^^"-^)^/^™ + . . . 



(17) 



where . . . denotes terms subleading to the ones showna It is straightforward to show that 

7;, = 2(A$-4), 7^ = min{2(A$-4),2(A,-4)} . (18) 

Also, Ba and bh can be determined algebraically in terms of a,,, 0$, A^,, and A4,. Using (flTI) one 
can express v in terms of bh, 

v-vin = \bue^^^* , (19) 
to leading order in t> — vib,. Also at leading order, one has for v ^ f ir 

^drag . (20) 

Comparing (fT9l) and (l20l) . one arrives at (lH). 



B Some details on obtaining the expression for ^j^j^. 

To prove ([5]) consider a positive but low temperature background as discussed in the main text, but 
with a vector perturbation of the form 

A.. = a(r) g,,. = e2^M^?(r) . (21) 

The equation of motion for g and a is a set of two coupled second order differential equations. As 
discussed in [|20|, one of the solutions to this set of equations is a boost of the background metric 

a = —(pu g = {h — l)u (22) 

with u an infinitesimal boost parameter in the direction (we use x = (a;^ , , x^)). The explicit 
form for the equation of motion for g can be written as a total derivative by making use of the 



"^There can also be an r-independent constant term in the expansion of A which can be removed by gauge fixing ; 
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equations of motion for the background. After integrating it once it takes the form 

e^Ag, ^ ^2A^^, ^ (23) 

where Q2 is an integration constant. Expanding 

a = -iiv, + \{J.)e-^^ + ... $ = + \{Jt)^~^^ + ■■■ 9 = -\{Tt.)e-^^ + • • • (24) 

with (J^) and (T^,/) the boundary theory charged current and stress tensor, we can obtain Q2 
in terms of boundary theory quantities. Using the relativistic two-fluid model for superfluids, 
introduced in ll22l and elaborated upon in [SI |23l |20l, it is possible to express (J^j) and (T^,/) in 
terms of the hydrodynamic quantities: s, T, /i, p^^, the velocity of the normal component u and 
the superfluid velocity v. This leads to Q2 = sTu.j.i + (n^i — Vxi)ppn- For the particular boosted 
solution (I22I) where u^i = v^i = u, (|23] ) reduces to 

e^^h' - e^H^' = (25) 

with Qi = sT. The identity (l25l) can also be obtained through a Noether charge argumentjf] Using 
(l23l) and (|25|) the equation of motion for a may also be integrated: 

he^^ia^' - a'4>) = hQ2 - gQi + Q3 (26) 

where Q3 = (vx — Ux)iJ.pn is another integration constant which has been evaluated by expanding 
(|26|) near the AdSs boundary. 

Next we evaluate (|26l ) in the deep JR. There, 

= 0' = O, h = v^^, g = {v^^-l)vx. (27) 

The first two expressions follow from (fTTI ) and the last expression follows from solving (|26l ) and 
(|23] ) with Qi = Q2 = Qs = and realizing that at zero temperature there is no normal component 
so that the only possible flow is given by the boosted solution (|22l) with the velocity of the super- 
fluid component being the boost parameter. As discussed in the main text, at low but non- vanishing 
temperatures one finds that (p, h and g approach the values (|27]) in the deep IR, before h vanishes 
at the black hole horizon. Thus, evaluating (|26l) for a positive but low temperature solution in the 
region where (ITTI) is valid, we obtain ([5]). 

^We thank A. Nellore for explaining the Noether charge argument to us. 
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C The stochastic force 



To study the stochastic forces and how they lead to ([8]), we consider perturbations of the string 
around the trailing string shape: 

/ r + C(r) \ 

VT + vC,{r) + ^{r) + 5x^{t, r) 
dx^ij, r) 
5x^{t, r) 
r 



(28) 



Plugging this ansatz into the action Q and expanding to quadratic order in the 5x^ leads to 



-"String 



2Tia'\fh 



e2Ah2^,2 _ ^{da5x^f - ^ ^{da^x'f + ^(^X^) , (29) 



i=2,3 



up to total derivative terms, where by {daSxY we mean 7^^* daSx dbSx, and 



Kr 



.fiA-2A, 



2Tia' h 



K' 



T 



h Q'^, 



2-Ka' ^/hlQ^ 



(30) 



Note that Kl and Kt are both finite and smooth at the worldsheet horizon. In what follows we 
would like to treat the transverse and longitudinal directions simultaneously, so we will use a 
subscript M = L,T where needed. The directions and x^ are associated with M = T while the 
x^ direction is associated with M = L. 

Because neither Km nor 7"'' depends on r, we may express the general solution to the linearized 
equations of motion following from (|29l ) as 



(31) 



6xHT.r)= I 0^(cu)e-^'^>M(^,r) 



where by assumption tjjM{uj, r) ^ 1 as r ^ 00 (the boundary of AdS5). At the string worldsheet 
we require that tjjM has only infalling modes. The equation satisfied by i>Mi^, r) is 



(32) 



This differential equation can be solved approximately for small uj using the method of matched 
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asymptotic expansions. Briefly, the method is to match the far-field solution 

r) = 1 + CU^) / ^ , (33) 

which solves (|32l) to order 0{uj'^), with the near-horizon expansion 

^M'^'i^, r) = Cnear(^)(r - r,)— /^^^^^ . (34) 

The integration constants Cfai and Cncar are r-independent, and in (|34|) we have chosen the infalling 
solution, associated with the retarded two-point function. Matching the near horizon behavior of 
(1331 ) to the small cu expansion of (|34l) leads to 

CU^) = -4^ {Km drYl l-r + ^(^') • (35) 

Given a solution ipM{^,r) to (|32|) . there is an established way to extract both the two-point 
functions of the stochastic forces on the probe and the strength km of these forces in the low- 
frequency limit [|25ll26ll27ll28l . The symmetrized Wightman two-point function is 



Gm{uj) = - coth ( -— lmG\'}iiu) = coth -— ) lim Imi^A^.i^MT'^'^r.^Af} • (36) 



r— >oo 



Here G^^}{uj) is the retarded two-point function. Plugging ( |33] ) into (|36l) gives 



Gm{uj) = - coth ( -— ] Im CiU^) + 0{uj) . (37) 



The strength of the stochastic forces is now given by 



= lim Gm{u;) = {KMdrYl L_, • (38) 



Taking the u ^ limit means that we are approximating G'j\/(t) by KMS{t). Obviously, this 
is valid only in the limit where we integrate the Langevin equation over a time longer than the 
characteristic timescale of Guit), which is 1/Tws- 

Using the expression for the worldsheet metric discussed in the main text and equations (l30l) 
dl]) and ©, one obtains the expression for kt in ([8]). From (|30l) and (l38l) one can see immediately 
that 

«L = 



Kt 



= ^,1:^ • (39) 
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To simplify this we start with (fT2l) and note that 
Which leads to in ([8]). 
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